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1. Introduction 



In the study of reaction-diffusion processes non-hermitian chains appear in a natural way. 
However, their properties have not yet been studied extensively. For instance, the effect of 
boundary conditions on the spectrum of non-hermitian Hamiltonians is unknown, although a 
rich and surprising structure arises here. In this article we study an asymmetric XXZ-chain 
with modified boundary terms and we show that these boundary terms give rise to new and 
interesting features. 

In this article, two examples are treated that both show the appearance of boundary- 
induced phase transitions. As a first example, we study the Hamiltonian of the asymmetric 
diffusion model on a one-dimensional lattice with open boundaries, injection and ejection of 
particles at the edges of the chain. For this model, all known results are valid only for the 
stationary state. The time evolution is studied for the first time in this article. In this model, 
particles of one species can diffuse on the lattice with length N with a rate p to the right and 
q to the left. Particles are injected at the left boundary with a rate a and extracted with a 
rate /3 at the right boundary. The Hamiltonian of this model is non-hermitian and is given 
by the following expression: 

N-l 



3=1 



■B 1 +B N (1.1) 



with 



a 

Bi = -(^-2^ + 1) 

Bn = -§K + 2<r+-l) (1.2) 

The bulk term of this Hamiltonian corresponds to the asymmetric XXZ-Hamiltonian with 
anisotropy A = This model and versions of it have been studied extensively. The total 
asymmetric diffusion model (with q = 0) was first defined by Derrida, Domany and Mukamel 
]l| who found the phase diagram for the current in the stationary state as a function of a and f3 
by application of mean-field theory. The ground state of the model can be determined exactly 
by a matrix product formulation ||, |J. Two massive and one massless phase have been found 
for different regions of the a- (3 plane. The massive phases are separated by a line on which 
the spectrum is massless. In this context massless means that the spatial correlation length 
diverges, which implies an algebraic falloff of the spatial concentration profile. On the other 
hand, massive means that the spatial correlation length remains finite. The phase diagram 
for the partially asymmetric system (p ^ q ^ 0) was presented for the first time in Essler 
and Rittenberg || calculated correlation functions for this system with the help of the matrix 
product formulation. 

The case of periodic boundary conditions has previously been studied by Gwa and Spohn 
H for q = and p = 1 using Bethe-Ansatz techniques; for arbitrary values of p and q, results 



can be found in 0, §. 

The first part of this article deals with the analysis of the spectrum of H for q = 0. 
The Hamiltonian is known to be integrable |J, but because of the lack of a reference state, 
the model cannot be treated by a Bethe-Ansatz. Therefore numerical methods are applied 
here. We were interested in determining the correspondence between the spatial correlation 
lengths already known from previous work described above and the dynamical properties 
of the Hamiltonian, i. e. the time correlation length r given by the inverse of the smallest 
energy gap Eq with respect to the ground state of the Hamiltonian [TO . Eq is determined 
for lattice lengths of 2 < N < 18 sites and then extrapolated to the thermodynamical limit. 
We find that r remains finite in the massive phase and diverges algebraically in the massless 
regime. The boundary induced phase diagram known for the stationary state is reproduced 
by the spectrum of the Hamiltonian; so the dynamics of the system exhibits the same physical 
properties as the stationary state. 

The second part of the paper is devoted to the antiferromagnetic version of the 
Hamiltonian —H where H was mentioned above. It can no longer be interpreted as 
Hamiltonian of a reaction-diffusion sytem, but nevertheless it describes interesting physical 
models, e.g. imposing periodic boundary conditions it can be viewed as the logarithmic 
derivative of the transfer matrix of the six-vertex model JTT[ and it can be used to describe 
the equilibrium shape of a crystal |12| or a surface growth model |13|]. The phase diagram of 
the antiferromagnetic chain with periodic boundary conditions has been determined recently 
14|j . A natural question arises whether different boundary terms alter the phase diagram in 



such a way that boundary induced phase transitions occur and indeed this phenomenon can 
be found again. It manifests itself clearly in the study of the partition function of the system. 
Before presenting our examples we give the definition of the partition functions we are going 
to use. For periodical systems the partition function is given by 

Z = lim tiz^ {H - eaaN) (1.3) 
for open systems by 

Z = lim tiz^ {H - foo - eocN) (1.4) 

Here £ is a normalization constant, the bulk free energy and /oo the surface free energy. 
For a comparison with the results for the asymmetric XXZ chain let us first repeat some 
known results for the symmetric XXZ-chain with different (but not lattice-length dependent) 
boundary conditions. It has been studied extensively by using numerical and analytical 
methods (Bethe-Ansatz) [|15|, I6[ . In the regime — 1 < A < 1 the spectrum is massless and 



can be described by a representation of the Virasoro algebra with central charge c = 1. The 
partition function corresponds exactly to the one of a Coulomb gas with electric charges n 
and magnetic charges m |T7[ where n corresponds to the total spin of the chain. It can also 



be described by a Gaussian field theory with compactification radius R [fL8||. It is given by 
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the following expression (using an appropriate normalization £(A): 

Z = z~& J2 z^ +R2m2) K(z) (1-5) 

mdZ,n 

with 

oo 

M^na-**)- 1 • (i- 6 ) 

8=1 

In eq. (|1.5| ) n takes integer values for an even number of sites and half-integer values for an 
odd number of sites. 

The spectrum of the chain without any boundary terms is given by a Coulomb gas with 



only electric charges n. Here the partition function reads [15 



1 n 2 

Z = z^U v {z) (1.7) 

nef 

This expression is valid for the open symmetric XXZ-chain as well as for the open asymmetric 
XXZ-chain since there exists a mapping of the symmetric onto the asymmetric XXZ-chain. 

For the periodic asymmetrix XXZ-chain that is normally massive in the antiferromagnetic 
region, a massless phase was found for A < 1 || |TJJ]. Here the operator content contains a 
term which is proportional to the lattice length N so that the partition function is the one 
of a modified Coulomb gas. It has been calculated by Noh and Kim using Bethe-Ansatz 
calculations M. 



Z = z~& £ z ^ +R2m2+2imNy) Ii 2 v {z) (1.8) 

mdZ,n 

The compactification radius R depends on A and q/p, n takes integer values for even lattice 
lengths and half-integer values for odd lattice lengths. 

These three examples show that boundary conditions can change the spectrum of the 
asymmetric XXZ-chain drastically. 

The boundary conditions we analyse in this article (with additional non-diagonal 
boundary terms) have not been treated before. Our new results are obtained by determining 
the lowest-lying eigenvalues of the Hamiltonian up to 18 sites in the general case and 21 sites 
in the CP-invariant case (a — f3). As the Hamiltonian does not have any symmetries in the 
general case and is non-hermitian, the diagonalization requires a large amount of CPU-time. 
A version of the deflated Arnoldi algorithm has been applied to quantum spin chains for the 
first time here. The Arnoldi algorithm that reduces to the Lanczos algorithm in the case of 
hermitian matrices was already used for the determination of the spectrum of the Potts model 



19 1 . Another example where a non-hermitian chain was treated numerically can be found in 



The analysis of the spectrum for p > q and a, (3 > suggests surprisingly the partition 
function of a Coulomb gas that has only magnetic charges. Additionally, the magnetic 
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quantum number m is shifted by an imaginary amount ix. 

Z = z -uJ2 z R2{ - m+l{x+Ny))2 H v {z) (1.9) 



mef 



Here x = x(|, ^, ^) and y = y(^). For even lattice lengths m takes half-integer values, for 
odd lattice lengths integer values. The spectrum is given by non-unitary representations of 
the U(l) Kac-Moody algebra. For q = the parameters y and R correspond exactly to the 
ones that appear for the periodic system ||. The term proportional to x however is caused 
by the boundary terms alone. For a(3 = the spectrum is massive, so we can conclude that 
the boundary terms give rise to the conformal invariant structure of the spectrum. 

We also studied a toy model consisting of a simplified version of the Hamiltonian (IT ) that 



can be diagonalized in terms of free fermions. The properties of the spectra that have been 
found numerically can be found again in the analytical results for the toy model. Additionally, 
we gain some insight into the role the different boundary terms play for the spectrum of the 
Hamiltonian. 

The article is organized as follows: In the first part, we define the total asymmetric 
diffusion model with open boundaries and additional injection and extraction terms. We 
present our numerical results for the temporal correlation length in the second part and 
compare them to the phase diagram of the stationary state and the known expressions for 
spatial correlation lengths. The third part deals with the antiferromagnetic version of the 
Hamiltonian. From the finite-size scaling behaviour of the lowest-lying energy levels we 
conclude the form of the partition function (|1.9|). The fourth part describes analytical results 
obtained by diagonalization of the toy model for non-hermitian boundary conditions. Here the 
characteristic properties of the Hamiltonian that has been treated numerically are reproduced, 
thus the calculations for the simplified model give similar results to our numerical analysis 
of the more complicated model. We close with a discussion of our results. In the appendix 
tables of the numerical results of the second and third part can be found. 



2. Asymmetric diffusion model with boundary terms 

We consider a model defined on a chain with N sites. Each site can be either occupied by a 
particle of species A or empty. For the dynamics, we consider only processes that involve two 
neighbouring sites. The following processes are allowed: 

diffusion to the right: A + — > + A with rate p (2 1) 

diffusion to the left: + A -> A + with rate q 
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Additionally, we allow processes that occur only at the two edges of the chain: 

.. . , , r , + 0^ A + ) . . 

miection at the left edge „ , > with rate a 

J to + A + A I 



A + A — > A -\- I 
extraction at the right edge: , > with rate 3 



(2.2) 



To each site j of the chain we attach a variable jj that takes the value 1 if the site is 
occupied by a particle and if the site is empty. A configuration of the model is described by 
a set of variables {7} = {71, 72, ... , 7tv}. The probability that a state (/x, v) on two adjacent 
sites will change into the state (p, a) after one unit of time is denoted by: 

r£; fr,u)?(p,a). (2.3) 

All reactions changing the state (/x, v) into any other state are summarized in the rate 

r^= £'r^\ (2.4) 

P,cr=0 

where the prime is always used to indicate that in the sum the case (p, a) = (/x, v) is excluded. 
This definition ensures the conservation of probabilities. 

The dynamics of the system is given by the time evolution of the probability P({ , ~f},t) 
to find the system in the configuration {7} at time t. The time evolution is described by the 
master equation [2~I| 

-§P({ 7 };t) = x![- r 7 fc ,7 fe+1 ^(7i,---,7L;t) 
01 k=i L 

+ E 'r|;^P(7i,...,7fc-i,^,^ + i,7fc+2,...,7L^)] • (2.5) 

<5fc,<5fc + i=0 

The master equation can be mapped on a Euclidean Schrodinger equation. For this 
purpose, one defines a ket vector describing the probability distribution in a 2 N dimensional 
vector space C® N : 

\P(t)>=J2P({l},t)\{l}> • (2-6) 

{7} 

The master equation then takes the form 

d_ 

Of 



P{t) >= -H\P(t) > . (2.7) 



where the Hamiltonian is a sum of nearest neighbor terms 

N-l 

H = y j H^k+i (2.£ 
fc=i 



The 4x4 matrices H k:k+1 are given by the following expression with respect to the basis of 
states 1 11 >, 1 10 >, |01 >, |00 > 



Hk. 



fc+i 



/ \ 

q -p 

-q p 

V o o o o j 



for k = 2, 3, ...,7V -2 



(2.9) 



The boundary terms which arise due to the processes of injection on the first lattice site 
and extraction on the last site can be written in this formalism as 



Hi, 
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(2.10) 



The Hamiltonian can also be expressed in terms of Pauli spin matrices a x , a v and o z : 

L-l 



H 



-E 



qa j a? +1 +paj~a j+1 + 



p + q 



1) + ^— (a: 



°3. 



+ J B 1 + J B JV (2.11) 



with boundary terms 



a 



Bi= + 
B N = -^(a z N + 2a+-l) 



(2.12) 



The spectrum of H stays invariant if a and /3 are permuted. This permutation can 
be obtained on the level of the Hamiltonian by applying the similarity transformation 
a ± — > cr T ,a z — > — <r z and a simultanous reflection of the system. Physically one can 
understand this property by denoting that particle-hole symmetry holds: The model can be 
seen as describing particles moving in one direction or holes moving in the opposite direction. 
If a — (3, the Hamiltonian is invariant under this transformation (CP-symmetry). 

The Hamiltonian can also be mapped by a unitary transformation U onto the symmetric 
XXZ-chain with different boundary terms. This correspondence will prove very helpful in the 
sequel and therefore it will be given explicitly here. If we take 



U -- 

with Q 
H' 



N 
3=1 



u, = h 



1 o 

AO*" 1 



in 



(2.13) 



and arbitrary A the transformed Hamiltonian H' = U HU is: 

p j 



N-l 

E 



a 3 a 3+l 



V V 
a 3 a 3+l 



Q + Q- 1 



a 3 a 3+l 
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+ B[ + B' N . (2.14) 

This is the UqSU(2) invariant Hamiltonian |22|] with additional length-dependent boundary 
terms: 

B[= fK 2 -2Q^V + l) , 

B' N = -^{v z N + 2Crfo+-l) . (2.15) 

Although this Hamiltonian is known to be integrable ||, it cannot be solved by the Bethe- 
Ansatz because one cannot construct a reference state. The non-hermitian boundary terms 
induce completely different properties of the Hamiltonian H in comparison to H' without B[ 
and B' 2 . The reason for this can be understood by considering the matrix form of H. As the 
bulk terms (without B[ and B' 2 ) commute with the total spin S z = X/JLi °~* ! > these bulk terms 
can be written as a matrix in block-diagonal form, where each block operates in a sector with 
fixed total spin. The boundary terms change the total spin by ±1, so they will appear in 
blocks below or above the diagonal. Only in the case where a/0 and (3 ^ 0, they will give 
contributions to the spectrum. In all other cases the spectrum of H is the massive spectrum 
of the L^gS'?7(2)-invariant XXZ-chain. 



+ 



Q-Q- 



1 o~ 



3. Results for the total asymmetric diffusion model with boundary terms 

3.1. Analytical results for the stationary state 

For future reference, let us first present some of the results achieved previously for the total 
asymmetric diffusion model. Up to now, all results are known for the stationary state. In 
|| the phase diagram for the current and the spatial profile of the concentration have been 
determined as functions of a and (3. The current through the bond j for the configuration 7 
is defined as: 

Jk = J2ik(i-ik + i)P(h},t) . (3.1) 

{7} 

The phase diagram is given in figure 1. There are three main phases, characterized by 
the density of particles: The high-density phase A, the low-density phase B and the maximal 
current phase C. Phases A and B are further subdivided into A ly A 2 and B\,B 2 . Phases A 
and B are separated by a line which is called 'coexistence line'. In the thermodynamic limit, 
the currents in the three phases are given by: 

A : j = /?(!-/?) 



B : j = a{l — a) 



S 



(3.2) 



The density profile of the concentration on the stationary state obtained in [[J allows to read 
off the spatial correlation length £ in the different phases. In the pase A, £ is defined by 

< rik >= const. + const. exp(— k/£) ; (3.3) 
in the phase B it is given by 

.k-N-1. , , 

< rife >= const. + const. exp( ) . (3.4) 

In the different phases the expressions for the correlation length are: 

A 2 :? = -ln(4/3(l-/?)) 

= []»«)] - 
L a(l — a) J 

B 2 :£ = -ln(4a(l-a)) (3.5) 

In the phase C and on the coexistence line the one-point function < > shows an algebraic 
behaviour: 

1 l-Spi _ x 
C :< nAr_ fe > = 2 

2 2^ 

, 1 - 2a , . 

coexistence line :< n k > = a + k — — — (3.6) 

3.2. Numerical results for the time correlation length 

In this chapter we compare the phase structure of the stationary state to the spectrum of 
the Hamiltonian. We investigate the gap between the ground state and the first excited state 
of the spin chain. This energy gap Eq allows us to read off the time correlation length r 
directly: 

t-'-Eq . (3.7) 

We determine the lowest lying excitation energies for lattice lengths 2 < iV < 18 by 
diagonalizing the Hamiltonian numerically, using a version of the deflated Arnoldi algorithm 



23| that reduces to the Lancos algorithm in the case of hermitian matrices. The eigenvalues 
were then extrapolated with the help of the BST-algorithm [Q. a and /3 have been varied in 
steps of 0.1 between 0.1 and 1. The following form of the energy gap has been found in the 
different regions of the phase diagram. In the phases A\ and Bi, Eq is a function of a and 
(3, while in phase B 2 , Eq depends only on a and in phase A 2 only on (5. In phase C and on 
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the coexistence line the system is massless: 



a 



A 1 
A 2 
Bi 
B 2 
C 



E G ~- 

E G 
E G ~- 

E G = 

Eq ~ 
Eq ~ 



m(a, j3) 
-- m(J3) 
m(a, (3) 
- m(a) 
N~i 

N~ 2 



Here m denotes the mass of the spectrum. 

Table I shows extrapolants for the energy gaps. One sees clearly that in the phase B\ the 
energy gap depends on both a and (3 while it depends only on a in the phase B 2 . We did 
not give any extrapolants for the phases A\ and A 2 since the spectrum is symmetric under 
permutation of a and (3 as described in section 2. Table I also shows that the mass gap 
vanishes in the phase C and on the coexistence line. In Table II extrapolants for the exponent 
| in the phase C and for the exponent 2 on the coexistence line can be found. 

Comparing these results with (|3.5| - |3.6| ), we can see that the temporal correlation length 
shows the same behaviour as the spatial correlation length in the different phases: in the 
phase C and on the coexistence line r diverges algebraically while it remains finite in the 
phases A and B. The exponent | has also been found for the total asymmetric diffusion 
model with periodic boundary conditions ||. The exponent 2 was found for the symmetric 
XXZ chain with periodic boundary conditions at the point A = 1. The periodic model can 
be mapped onto a model for surface growth || [7], This mapping can be formulated 
analogously for the open chain with additional boundary terms that is treated here. Then, 
in the language of growth models, the exponent | describes KPZ-type growth [^6| while 
the exponent 2 corresponds to Edwards- Wilkinson growth behaviour [Tj|. In phases A 2 and 
B 2 t depends only on (3 or a respectively. However, the numerical values for the temporal 
correlation length are different from the ones for the spatial correlation length. The fact that 
the temporal evolution of the system reflects itself in the stationary state (for t — > oo) is not 
yet understood. This is not true in general, as there are examples for a different behaviour 
of the spatial and the temporal correlation length ||27|| . However, in this case it is a striking 
and unexpected feature that shows how strong boundary conditions change the properties of 
non-hermit ian quantum chains. 



4. The antiferromagnetic chain 

In the following, we will treat the antiferromagnetic Hamiltonian which can be obtained from 
the ferromagnetic one (that has been treated before) by applying a similarity transformation. 
As has been shown in chapter 2 the ferromagnetic chain can be connected to a XXZ-chain 



10 



with anisotropy A = > 1. The antiferromagnetic chain corresponds to A < 1 and can be 
written as 

J '~' - -!- - - I 1 '"",, - - 7 
tin rr ■ , -+- nci a , in n , — I ; — — ^ 



with boundary terms 

s? = -fte + 2*r + i) . 

^ = ^(^ + (-1)^-1) . (4.1) 

Due to the minus sign in the transformation rule, the low lying excitations of the 
ferromagnetic chain correspond to the highest states of the antiferromagnetic chain and vice 
versa. 

The numerical studies concentrate on the case q = and p = 1 with a, (3 > 0. For this 
choice of parameters the extrapolated values converge very well. 

The analysis of the spectrum reveals that it is conformal invariant. We studied the 
finite-size scaling behaviour of the model in order to determine the operator content of 
the underlying conformal field theory. For conformal invariant systems, the ground state 
is supposed to take the following form for finite lattices and for periodic boundary conditions 



28|, |29| 

E (N) ttCc 
N 6o ° 6iV 2 



o(N~ 2 ) (4.2) 



The excited state E r satisfies 30 



S r = lim £-{E r {N) - E (N)) = (A + r) + (A + f) ; (4.3) 

N-^oo 2711; 

for open boundary conditions, one gets an additional surface term in the ground state 
energy [j3l"H : 

Eo(N) /«, ixic 2 

The energy gaps scale here as 

S r = lim ^(E r (N) - E (N)) = (A + r) ; (4.5) 

N— >oo 7Tc 

Here ^ is a normalization constant, c the central charge of the Virasoro algebra. 

In our case, the spectrum had to be treated separately for even and odd lattice lengths. 
Since the ground state was calculated for an odd number of sites, we used interpolated values 
from odd lattice lengths for the values of the ground state for even lattice lengths. While 
the ground state energy is real for all lattice lengths, most of the excited states have a 
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non- vanishing imaginary part. This is a feature that did already appear for the periodic 
asymmetric XXZ chain H and in the calculation of the operator content of the five vertex 



model defined on an anisotropic lattice |32 



We considered the imaginary part of the energy gap 

J = lim ^-lm(E(N)) , (4.6) 

N^oo Tit 

the correction to the real part 
N 

Re(£) = lim —Re(E(N) - E (N))) (4.7) 

N^oo Jit 

and the correction to the imaginary part 
N 

lm(8)= lim —lm(E(N)-T) (4.8) 

N~*oo fit 

The data reveals the following finite-size scaling behaviour for the normalized eigenvalues: 
^-(E™-E' Q ) = R 2 (m + i(x + yN)) 2 + r , reN , (4.9) 
and the ground state behaves as 

E' = E + Tia^r(x + yN) 2 . (4.10) 

Here m is a quantized number and takes for odd lattice lengths the integer values m = 
0, ±1, ±2, ■ ■ •, for even lattice lengths half-integer values m = ±|,±|,.... The precise relation 
between the extrapolants and the constants R, x and y is given by: 

= 2R 2 my, 

Re(£ r m ) = R 2 m 2 + r, (4.11) 
Im(£ r m ) = 2R 2 mx. 

Data for r = is shown in table III. As H is real, we find for each complex eigenvalue E also 
the complex conjugate E*. Therefore all tables only show data for m > 0. Data for r > 
can be found in table IV for m = and m = \. The corresponding data for m = 1, |, 2, | 
and 3 has also been obtained. 

The next constant that has to be determined is the conformal charge d . The ground state 
energy in the thermodynamic limit, e^, is already known from the periodical system. The 
surface energy can be obtained by extrapolation of 

U = lim (E (N) - N eoo ) . (4.12) 

N— >oo 

The results are given in table V. Now c can be obtained via 
24iV 

c = lim -—(E (N) - N eoo - U . (4.13) 

N^oo Tit 
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The numerical values for c are not constant for different values of a and (3 (see table VI). 
However, we can shift all energy gaps, independently of the sector by a constant term 
24R 2 x(a, (3) 2 that depends on a and (3 (table VII). Absorbing this shift into the ground 
state, one can define a new central charge d 

d = c-2AR 2 x 2 . (4.14) 

which is indeed constant within the numerical errors. The numerical values show an excellent 
agreement with d — 1 (table VI). 

The term 7r£7™, the imaginary part of the energy gaps, takes (independently of a and (3) 
for odd lattice lengths multiples of the same constant 1.8854.... that already appeared in the 
calculation of Gwa and Spohn || for the periodic system as imaginary part of the smallest 
energy gap. Using Bethe Ansatz calculations for the first and second smallest eigenvalue in 
the sector with spin they obtained the following result 

eoo = 0.690140115. (4.15) 

(E 1 - E ) per = 6.5776787iV~ 1 + il.885456427. (4.16) 

This result can be used to obtain estimates for the parameters R and y, if one assumes 
that also the spectrum of the periodic chain is characterized by these parameters, and the 
finite-size scaling of the lowest lying state with spin is given by || 

N ,„ „ , 1 , „ 



} (Ex - E ) per = -R 2 + iR yN . (4.17) 
Zttc, 2 

The estimates that have been obtained using this assumption can now be compared to 
the numerical results. Tables III and IV show a comparison between the values obtained from 
the analytical calculation and the numerical results. The normalization constant £ has been 
taken from numerical Bethe Ansatz calculations for the periodical system for up to 80 sites 
@: 

£ = 1.64784392694623 . (4.18) 

The data shows that indeed the periodic system is characterized by the same parameters as 
the open system with additional boundary terms. However, the full operator content of the 
two chains is different. This effect comes from the different boundary terms only. 

The degeneracies of the energy levels for r = const, are described by the characterfunction 



of a U(l) Kac-Moody algebra This confirms the above result d — 1. The Kac-Moody 
algebra is defined by its commutation relations 

[T m ,T n ] = m5 m+n} Q , m,neZ . (4.19) 

The character function is given by 

X A, q (z, y) = tr(z Lo y T °) = z^U v {z)y' 1 . (4.20) 
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where L is a generator of the Virasoro algebra with conformal weight c = 1 that can be 
canonically obtained from the U(l) Kac-Moody algebra using the Sugawara construction 
||35|| . Here q is the eigenvalue of To, and the highest weight representation is [A, q\. A shift in 
the algebra characterized by a parameter (p 

T m = T m + (p5 mfi (4-21) 

does not change the commutation relations above but leads to a representation with highest 
weights 

A = -q 2 where q = q + (p (4.22) 

If if is chosen to be complex, one obtains a representation of the Kac-Moody algebra with 
negative conformal dimensions. In this case however T + ^ T and the representation is not 
unitary. 

In our case, the part of the energy corrections that is independent of N can be described 
by a non-unitary representation of a shifted Kac-Moody algebra. The parameter x depends 
on the boundary terms. The shift <p of the Kac-Moody algebra is given by <p = \[2%{x + Ny). 

The results for the antiferromagnetic chain can be summarized in the partition function 

Z = z -& z R2( - m+i ^ +Ny ^U v (z) (4.23) 

The operator content of this model corresponds to the one of a Coulomb gas with only 
magnetic charges and an additional term that depends on the lattice length N. All calculations 
for q ^ p for p > q show similar results but reveal that in this case x and y are functions 
of the ratio q/p. 

The case q > p has not been studied systematically. Here the spectra a purely real, but 
the convergence of the extrapolations is too bad to obtain precise estimates (the same holds 
for the case p = q, i.e. the symmetric XX Z chain with additional boundaries). 



5. Toy model 

In this chapter we study analytically a modified version of the Hamiltonian treated before. The 
results for different boundary conditions show that the main new features in comparison to 
the hermitian chain can already be found in this toy model. This model allows to understand 
which term in the Hamiltonian is responsible for the new and unexpected contributions to 
the spectrum. 

We will concentrate on 

2V-1 

For a = or (3 = the spectrum consists only of the N times degenerate eigenvalue zero. So 
the whole structure of the spectrum is caused by both boundary terms together. 
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This Hamiltonian can be diagonalized in terms of free fermions. Moreover H(a,/3) can 
be transformed into H(—a, —(3) by applying the transformation — > —a ± . For a[3 > the 
characteristic properties of the spectra described previously are reproduced. 

In order to write H in terms of free fermions, we have to obtain a bilinear expression in 
cx-matrices so that standard fermionisation techniques can be applied Technically, this 
can be achieved by appending one lattice site at each end of the chain, site and site N + 1 
PTfl . The Hamiltonian then reads 

N-l 

H'=J2 °t°i+i + aa o a ± + ^n^n+i (5-2) 
j=i 

As erg and <J% +1 commute with H' thus being constants of motion, the spectrum of H' 

decomposes into four sectors (++, H — , — h, ) corresponding to the eigenvalues ±1 of (Jq 

and <jfj + i- The eigenvectors of the extended Hamiltonain include those of H. Therefore we 
can obtain the eigenvectors of the original problem by projecting onto the (++)-sector. 
Defining new operators Tj~ and r~ |38| by 

= (IKM. ^ r = (IKK (j = 0,...,iV + l) (5.3) 

i<j i<j 

that obey the anticommutation relations of a Clifford-algebra 

{rf , rj} = 28% = 0, .., N + 1; = ±1). (5.4) 

one can rewrite H' as a bilinear expression in rf and tJ 

N-l 

H' = ~ E E A T^N+i + B^t$t v x + B^r N r v N+1 (5.5) 

H,i/=±1 j=l 

with 

Here we chose as a basis 

' - ( aZ a a Z ) ■ (") 

Then a second linear transformation 

N 

= E E MOJr" (7 = ±1) (5.8) 

3=0 jt=±l 

with 

{If, IJ} = (i, j = 0, .., JV + 1; //, v = ±1). (5.9) 

yields 

N+l 

ff' = E An ^T"T+ (5.10) 



n=0 
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The eigenvalues A n and the vectors ^ = ((ipZ)o , (^n)o , Mjv+i, (#»)jv+i) are g iven b Y 
the solution of the eigenvalue problem 



M(V>+±#-) = T A n (^+±^-) . 

where M is a (2N + 4) x (2N + 4) matrix which is given by 

/ B \ 
-B T A 

-A T A 



(5.11) 



M 



\ 



—A T C 



(5.12) 



J 



The digonalization of M leads to one eigenvalue and the following equation determining 
the other eigenvalues: 

(2A) 27V+2 = 4a 2 (3 2 (-l) N . (5.13) 
leading to the solutions 

A = 0, 2A n = (2a/?)s£r(sin(fc) -icos(Jfe)) (5.14) 

with 

fe = 2n - 1 7T (0<n<N + l) . (5.15) 

Because the operators iT~T+ always have eigenvalues ±1 it is sufficient to determine the 
eigenvalues A n with a positive real part. The energy gaps of the spectrum are given by 2A n . 
For the calculation of the exact ground state, one has to sum all levels with negative energy 
(so that the fermi see is filled), which leads to 



N+l 

E = - J2 An 

n=l 



(2a/3) 



N+l 



2sin( 



(5.16) 



2N+2' 



5.1. Scaling of the energy levels 
For large N the ground state reads 

_ N l+hU>/3) , 7T 

— -C/n ~ 

7T 7T 24N 



' 121n 2 (2a/3)\ ... 
1 + ^— ^ + o(N~ 2 ) 



7T Z 



(5.17) 



For the study of the low-lying excitations one has to develop the expression ( |5.14j) around the 
minimum of the dispersion relation which is obtained for momentum ±|. For < n € iV 
one obtains: 



2A r . 



2n - 1 7T 
2 iV 



-i 1 



ln(2a/3)~ 
iV 



for 



7T 



(5.18) 
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2A jV+2 _ w « — ^— - + 1 f 1 + ^ j for fc>-. (5.19) 

Up to now we always treated the chain where two additional sites, and N + 1 have been 
added to the starting Hamiltonian. We obtain the spectrum of the Hamiltonian without these 
sites by projection onto the sector where the cx-matrices acting on the additional sites have 
eigenvalue 1. It can be shown that the spectrum of H with respect to the ground state Eq for 
large N is obtained for even N by combining only an odd number of fermions with energies 
2A„ and for odd lattice lengths by combining an even number of fermions. The ground state 
of the system is found to be in the sector with odd lattice lengths. 
For the normalization constant, one reads off directly £ = 1. 

We turn now to the determination of the partition function because this can be compared 
directly to previous results. Writing H = —(H — E ) and using the triple product identity 
|T8l we obtain 

tiz 6 = n£°=i(l + z n -h ia )(l + z n ~h' ia ) 

= E^ e | ^ m+ *t) 2+ 4 rc=i(i - z»y x (5.20) 

with 

a „ r ln(2aB) 1 , . 

-= x + y N, x = — — — , y = — . 5.21) 

2 y ' 2tt ' y 2tt V ; 

Here we obtain m e Z for odd lattice lengths and m G Z + ^ for even lattice lengths. 

2 

Absorbing the term ^- into the ground state energy, we get the following expression: 

This corresponds to a system with conformal charge c = 1. For the partition function we 
obtain 

oo oo 

Z = £ z 2{m+i{ X+ Ny)f (1 _ ^-1 (5,33) 

mef n=1 

This is exactly the same result we found for the more complicated model with R 2 = 2. 

Comparing the partition function for the toy model with the ones for the open and for 
the periodic XXZ-chain given in the introduction, we get some insight into the nature of the 
different new terms in the partition function. For asymmetric XXZ-chains without periodic 
boundary conditions, we get the spectrum of a Coulomb gas with one kind of charges only. In 
the case of the completely open chain we get electrical charges, in the case of non-diagonal, 
non-hermitian chains we get only magnetic charges. These boundary conditions also give rise 
to an imaginary shift ix in the magnetic charges that also appears for the periodic chain with a 
real twist |L5|]. In our case x depends on a, (3 and the ratio p/q. The lattice length-dependent 



term yN appears for all non-hermitian chains we mentioned in this article except for the 
completely open chain which is equivalent to the hermitian symmetric open XXZ-chain. 
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6. Conclusions 



In this article, the spectrum of the asymmetric XXZ-chain with non-hermitian boundary 
terms has been studied. Numerical methods have been applied for the determination of the 
eigenvalues. Analytically, we studied a toy model that reproduces the characteristic properties 
of the full Hamiltonian. In two different cases, we concluded that the boundary terms are 
responsible for the phase transitions of the Hamiltonian. 

We studied first the ferromagnetic chain with a, (3 > that describes the time evolution of 
a reaction-diffusion system with asymmetric diffusion in the bulk and additional injection and 
extraction terms at the boundaries. For the completely asymmetric chain the behaviour of the 
smallest energy gap that corresponds directly to the inverse temporal correlation length shows 
the same behaviour as the spatial correlation length in the stationary state determined by 
0. In phases A and B where the spatial correlation length stays finite, i.e. the concentration 
of particles shows an exponential decay in the spatial direction, we find a massive spectrum 
where the mass depends on the same parameters as the spatial correlation length in the 
different phases A\,A2, B\ and B2. In this way the subdivision of the phases A and B is also 
valid for the dynamical properties of H . However, the numerical values of the spatial and 
the temporal correlation lengths do not coincide. The phase C, where the spatial correlation 
length of the stationary state decays algebraically, exhibits a decay of the temporal correlation 
length with an exponent |. This exponent has already been identified for the model with 
periodic boundary conditions. On the coexistene line, we find a decay with an exponent 2. The 
deep reason why the phase diagram of the stationary state is reflected in the time-evolution 
of the system is not yet understood. 

The antiferromagentic chain for A < 1 shows a completely different behaviour. For a = 
or (3 = the spectrum corresponds to that of a XXZ-chain with additional cr z -terms at the 
boundaries and therefore is massive. Only when both a and j3 are non-zero, the spectrum is 
massless and can be described by a representation of a non- unitary U(l) Kac- Moody algebra. 
However, one term proportional to the lattice lenght arises in the finite-size scaling spectrum. 
The result can be summarized in the partition function of a modified Coulomb gas with only 
magnetic charges: 

Z = £ Z K 2 (™+i(*+Kv)) 2 n v (z) (6.1) 

mef 

The parameters y and R also appear for periodical boundary conditions while x is induced 
by the boundary terms. The analysis of the toy model that was diagonalized in terms of 
free fermions reproduces the same structure of the spectrum. The length-dependent term 



seems to be a common property of the anisotropy of the spin chain |32j and also appeared 
for the periodic chain || . The imaginary shift of the magnetic charge is an effect of the non 
hermitian and non- diagonal boundary terms. The influence of hermitian boundary terms 



and asymmetric interactions in the bulk will be the subject of a future publication fl39| . 
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The question arises now how the ordinary Coulomb gas model has to be modified so that 
imaginary and lattice-length dependent contributions to the operator content arise. Is it still 
possible to find a field theory that reproduces the spectrum of the modified Coulomb gas? A 
first step in this direction could be the calculation of the correlation functions that normally, 
for conformal invariant systems, exhibit critical exponents which are directly related to the 
eigenvalues of the corresponding Hamiltonian. Wether this relation is still valid in the case 
of non-hermitian chains still has to be clarified. 

The new and unexpectedly interesting structure that appeared in the two examples treated 
in this article suggests further research in the field of non-hermitian Hamiltonians. 
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Figure 1. Phase diagram for the total asymmetric diffusion model 
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Table I. Asymmetric diffusion model: Dependence of the extrapolated gap Eg on a and (3. Negative 



values indicate that errors are of the order of 10 



0\a 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.0 


0.12573 


0.04014 


0.01017 


0.00119 


0.00000 


0.00001 


-0.0000 2 


-0.00001 


-0.00000 


0.00000 


0.9 


0.12575 


0.04014 


0.01021 


0.00118 


-0.00003 


0.00000 


-0.00 004 


-0.00001 


-0.00000 




0.8 


0.12573 


0.04013 


0.01019 


0.00118 


-0.00001 


0.00000 


-0.00 004 


-0.00004 






0.7 


0.12587 


0.04012 


0.01019 


0.00114 


-0.00003 


-0.00002 


-0. 00004 








0.6 


0.12017 


0.04016 


0.01030 


0.00120 


-0.00001 


-0.00001 










0.5 


0.09999 


0.03331 


0.00861 


0.00110 


0.00000 












0.4 


0.07213 


0.02018 


0.00330 


-0.00005 














0.3 


0.04173 


0.00678 


-0.00001 
















0.2 


0.01428 


-0.00001 


















0.1 


-0.00000 
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Table II. Asymmetric diffusion model: Extrapolated exponents of the first excited state. 
f3\a 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



1.0 1.5000 1.470 1.456 1.480 1.4999 1.522 

0.9 1.499 1.452 1.503 1.4994 1.4999 

0.8 1.4997 1.430 1.436 1.499 

0.7 1.462 1.405 1.499 

0.6 1.434 1.377 

0.5 1.49999 

0.4 2.01 



0.3 1.98 
0.2 2.010 
0.1 1.992 



Table III. Extrapolated and predicted values of the energy gaps, computed from the finite size spectrum 



of 


-H for a = (3 


= 0.5, q = 0,p = 


1 and r = 0. 






in 


R 2 m 2 


2R 2 


my 


2R 2 mx 




prediction 


extrapolated 


prediction 


extrapolated 


extrapolated 


l 

2 


0.317648 


0.31763(2) 


0.182104 


0.182104 


-0.120 70(3) 


1 


1.270 594 


1.270 (6) 


0.364209 


0.3642(1) 


-0.2414(1) 


3 
2 


2.858836 


2.858 (9) 


0.546 313 


0.546 313 


-0.36210(7) 


2 


5.082375 


5.082 (3) 


0.728418 


0.728 (4) 


-0.4828(2) 


5 
2 


7.941211 


7.9(4) 


0.910522 


0.9(1) 


-0.60(4) 


3 


11.435343 


11.4(8) 


1.092 627 


1.0(8) 


-0.72(3) 
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Table IV. Data corresponding to (4.13) from the extrapolation of the spectra of —H for a = (i = 0.5, 
q = 0, p = 1 with 2 < N < 21; predicted degeneracies are given in square brackets; the predicted values 
of 2R 2 mx are determined from the extrapolation of the value 2R 2 mx for m = \ and r = 0. 







m = 








m= i 




r 


R 2 m 2 + r 


2R 2 mx 


2R 2 my 


r 


R 2 m 2 + r 


2R 2 mx 


2R 2 my 


1 [11 

L J 


1.000000 






[11 

L J 


0.31763(2) 


-0.120 70(3) 


0.182104 


2 [21 


1.999 99(7) 


0.0000(0) 


0.000 00(1) 


1 [11 


1.31765(7) 


-0. 120 702 


0.182104 


3 [31 

L J 


3.0000(0) 






2 [21 

L J 


2.317(4) 


-0.120 (7) 


0.1821(1) 




2.999 (8) 


-0.000(1) 


-0.000(0) 




2.3176(8) 


-0.120 70(2) 


0.182104 


4 [51 


4.0000(0) 






3 [31 

L^J 


3.3(2) 


-0.12(0) 


0.182(1) 




4.00(0) 


0.000(1) 


-0.000(0) 




3.317(7) 


-0.1207(0) 


0.182 1(0) 




4.0(0) 


0.00(3) 


-0.00(0) 




3.31(7) 


-0.1206(8) 


0.1821(1) 


5 [7] 


4.99(6) 






4 [5] 


4.3(2) 


-0.1(3) 


0.18(4) 




5.00(0) 


0.00(0) 


0.000 (0) 




4.31(8) 


-0.120 (7) 


0.182(1) 




5.0(0) 


0.0(0) 


-0.00(1) 




4.317(7) 


-0.120 (8) 


0.1820(8) 












4.31(6) 


-0.12(2) 


0.18(2) 


6 


5.99(7) 








4.31(7) 


-0.11(8) 


0.182 (0) 




6.0(1) 


0.0(0) 


-0.00(1) 












6.0(0) 


0.00(1) 


-0.00(0) 


5 [7] 


5.31(7) 


-0.121(4) 


0.1820(8) 




6.0(0) 


0.0(0) 


-0.00(1) 




5.3174(8) 


-0.120 (7) 


0.182(1) 












C O 1 ( A \ 

5.31(4J 


— 0.1z(zj 


0.18(2) 


7 


7.0(1) 








5.31(7) 


-0.12(1) 


0.18(2) 




6.99(9) 


-0.0(0) 


-0.000 (0) 




5.3(2) 


-0.1(3) 


0.18(4) 




7.0(1) 


0.0(0) 


-0.00(1) 




5.3(2) 


-0.11(5) 


0.1(8) 


8 


8.0(0) 






6 


6.3(1) 


-0.11(7) 


0.182(1) 




8.0(2) 








6.3(4) 


-0.1(1) 


0.18(1) 




8.0(1) 


0.0(0) 


-0.00(1) 












8.(1) 


0.(0) 


-0.00(5) 












8.0(1) 


0.0(0) 


-0.00(1) 










prediction 


0.000 000+r 


0.000 000 


0.000000 




0.317648+r 


[-0.120 70 ] 


0.1 82104 
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Table V. Dependence of the free surface energy on a, 0. 



0\a 


1.5 


1.3 


1.1 


0.9 


0.7 


0.5 


0.3 


0.1 


0.4959(9) 


0.30(9) 


0.123(1) 


-0.06(4) 


-0.25(6) 


-0.4631(0) 


-0.71(0) 


0.3 


0.9086(1) 


0.722(0) 


0.535(8) 


0.348(6) 


0.1564(6) 


-0.050216(5) 




0.5 


1.15624(0) 


0.96957(1) 


0.78344(0) 


0.596234(5) 


0.4040372 






0.7 


1.362(8) 


1.176(2) 


0.9900(9) 


0.8028(9) 








0.9 


1.555(0) 


1.368(4) 


1.18(2) 










1.1 


1.742(2) 


1.555(6) 












1.3 


1.928(3) 















Table VI. Extrapolated central charges. The value of R 2 x is determined from Eq 



a = 


loo 


R 2 x 


c c 






0.1 


-1.12(3) 


0.793(8) 


12.9(8) 1.0(8) 






0.2 


-0.6(2) 


0.506(0) 


5.8(3) 1.0(0) 






0.3 


-0.297(8) 


0.3339(8) 


3.1(1) 1.0(1) 






0.4 


-0.0336(9) 


0.2126(4) 


1.84(8) 0.99(4) 






0.5 


0.19737 


0.1207(0) 


1.2752(8) 1.0000(8) 






0.6 


0.40965 


0.0486(5) 


1.044(4) 0.999(7) 






0.7 


0.610(7) 


-0.008(7) 


1.00(0) 0.99(9) 






0.8 


0.8049(7) 


-0.054(6) 


1.04(8) 0.99(2) 






0.9 


1.1829(2) 


-0.0913(7) 


1.2(7) 1.1(2) 






1.0 


1.369(5) 


-0.1(2) 


1.3(8) 1.1(1) 






Table VII. Dependence of R- 


! x on a and 0. 






0\a 


1.5 


1.3 


1.1 0.9 0.7 


0.5 


0.3 


0.1 


-0.29(9) - 


-0.30(9) - 


0.32(5) -0.35(2) -0.39(3) 


-0.45900 


-0.565(6) 


0.3 


-0.06(8) - 


-0.07(8) - 


0.09(4) -0.12(1) -0.1(6) 


-0.22735 




0.5 


0.0379(7) 


0.02806 


0.01165 -0.01467 -0.05601 






0.7 


0.102(5) 


0.092(6) 


0.076(2) 0.049(9) 






0.9 


0.14(3) 


0.13(3) 


0.11(7) 






1.1 


0.17(0) 


0.16(0) 








1.3 


0.18(6) 
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